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Abstract 

In this article, we investigate the contribution of the high twist Feynman diagrams to the large-py 
pion production cross section in proton-proton collisions and we present the general formulae for 
the high and leading twist differential cross sections. The pion wave function where two non-trivial 
Gegenbauer coefficients ci2 and 04 have been extracted from the CLEO data, two other pion model 
wave functions, P2, P3, the asymptotic and the Chernyak-Zhitnitsky wave functions are used in the 
calculations. The results of all the calculations reveal that the high twist cross sections, the ratios 
R, r, the dependence transverse momentum pt and the rapidity y of pion in the <&cleo{ x ->Q 2 ) 
wave function case is very close to the Q asy {x) asymptotic wave function case. It is shown that the 
high twist contribution to the cross section depends on the choice of the meson wave functions. 
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I. INTRODUCTION 



During the last few years, a great deal of progress has been made in the investigation 
of the properties of hadronic wave functions [1-1 2]. The pion wave functions (also called 
distribution amplitudes -DA) [1] play a key role in the hard-scattering QCD processes be- 
cause they encapsulate the essential nonperturbative features of the pion's internal structure 
in terms of the parton's longitudinal momentum fractions Xj. Meson wave functions have 
been extensively studied by using QCD sum rules. The original suggestion by Chernyak 
and Zhitnitsky of a "double-humped" wave function of the pion at a low scale, far from the 
asymptotic form, was based on an extraction of the first few moments from a standard QCD 
sum rule approach[5], which has been criticized and revised in Refs. [6,7]. Subsequently, 
a number of authors have proposed and studied the modified versions of meson [7,8] and 
baryon wave functions [9,10]. Interesting nucleon functions have been constructed in the 
context of the quark-diquark model[ll]. Additional arguments in favour of a form of the 
pion wave functions close to the asymptotic one have come from the analysis of the tran- 
sition form factor 77* — > 7r° [12]. The measurements of this form factor by the CLEO 
collaboration are consistent with a near-asymptotic form of the wave function[13]. In [14], 
the leading-twist wave function of the pion at a low normalization point is calculated in 
the effective low-energy theory derived from the instanton vacuum. These results for the 
pion wave function at the low normalization point are close to the asymptotic form and 
consistent with the CLEO measurements. The authors have obtained a shape substantially 
different from the Chernyak- Zhitnitsky one because they have chosen a significantly smaller 
value of the second moment, and, more importantly, they have taken all the moments of 
the wave function into account. Their results support the conclusions reached previously in 
Refs. [6,7]. The QCD factorization theorems predict that the hadron-hadron cross section 
can be obtained by the convolution of parton distribution functions and a cross section of 
the corresponding hard scattering subprocess. The parton distributions are nonperturbative, 
process-independent quantities, which are specific to any given hadron. The hard scatter- 
ing cross sections are independent of all long distance effects and can be found by means 
of pQCD. In the framework of pQCD, the higher order corrections to the hard scattering, 
and therefore to the hadron-hadron process cross sections, have been calculated [15]. These 
corrections are large and change the leading order results considerably. Other corrections 
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to the hadron-hadron process cross sections and its different characteristics come from the 
higher twist (HT) terms. Thus, explicit HT effects associated with a meson bound state 
in the process nN — > 7*X and tiN — > 7X and their influence on the decay angular distri- 
bution have been found by Brodsky and Berger in Ref. [16]. It is important to note that 
the term "twist" ("twist"- means dimension minus Lorentz spin) is one of the characteristics 
of composite operators that occur in operator product expansion (OPE) in deep inelastic 
scattering (DIS). In DIS, the higher twist (twist >2) terms in OPE are associated with power 
suppressed corrections. On the other hand, in the various hard hadron-hadron scatterings- 
i.e. in the Drell-Yan process lacking an (OPE) description- the higher twist terms refer to 
contributions which are suppressed as 0(1/ Q 2 ), relative to the scale of the hard scatttering. 
Therefore, in order to calculate the power suppressed corrections to the hadron-hadron colli- 
sions cross sections, the Feynman diagram approach should be used. Indeed, in the context 
of this method, the factorization of the HT contributions of 0(1/ Q 2 ) in hadron collisions 
has been proven [15]. It is well known that in hadron-hadron scattering at the 0(1/Q 4 ) 
level, the factorization fails, which has been demonstrated by the existence of non-cancelling 
infrared divergences at two loops and by finite terms at one loop [17]. So, the extension of the 
factorization to 0(1/ Q 2 ) corrections in a large class of hadronic processes helps to place their 
existing treatments [16,18] on a solid foundation. In [15], the leading 1/Q 2 corrections to the 
Drel-Yan cross-section have also been obtained. It is worth noting that the normalization of 
the 0(1/Q 2 ) longitudinal structure functions in the Drell-Yan cross-section are determined 
by higher twist longitudinal structure functions in DIS. In other words, the 1/Q 2 corrections 
to the Drell-Yan process can be expressed in terms of the same multiparton correlations as 
in the DIS one. 

On the other hand, the results obtained are consistent with the lowest order calculation 
of the pion-hadron scattering carried out by Brodsky and Berger in Ref. [16]. By taking 
these points into account, it may be asserted that the analysis of the higher twist effects on 
the dependence of the pion wave function in pion production at proton-proton collisions are 
significant in both theoretical and experimental studies. 

Another important aspect of this study is the choice of the meson model wave functions. 
In this respect, the contribution of the high twist Feynman diagrams to a pion production 
cross section in proton-proton collisions has been computed by using various pion wave 
functions. Also, the leading and high twist contributions have been estimated and compared 
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to each other. Within this context, this paper is organized as follows: in section |Hj we 
provide some formulae for the calculation of the contribution of the high twist diagrams. In 
section llll| we provide the formulae for the calculation of the contribution of the leading 
twist diagrams and in section II VI we present the numerical results for the cross section 
and discuss the dependence of the cross section on the pion wave functions. We state our 
conclusions in section 

II. CONTRIBUTION OF THE HIGH TWIST DIAGRAMS 

The high twist Feynman diagrams, which describe the subprocess qi + q 2 — > 7i + (tt~) + 7 
for the pion production in the proton-proton collision are shown in Fig.l. In the high twist 
diagrams, the pion of a proton quark is directly observed. Their 1/Q 2 power suppression 
is caused by a hard gluon exchange between pion constituents. The amplitude for this 
subprocess can be found by means of the Brodsky-Lepage formula [19] 

M(s,i) = / dxi I dx 2 5(l - Xx - x 2 )§ 7T (x 1 ,x 2 ,Q 2 )T H (s,i;xi,x 2 ). (2.1) 
Jo Jo 

In Eq.(2.1), T# is the sum of the graphs contributing to the hard-scattering part of the 
subprocess. The hard-scattering part for the subprocess under consideration is q\ + q 2 — > 
(Q1Q2) + 7, where a quark and antiquark form a pseudoscalar, color-singlet state (qiq 2 ). 
The hard-scattering amplitude Th(s, t; x%, x 2 ) depends on a process and can be obtained 
in the framework of pQCD, whereas the wave function ^^(xi, x 2 , Q 2 ) describes all the non- 
perturbative and process-independent effects of hadronic binding. The hadron wave function 
gives the amplitude for finding partons (quarks, gluons) carrying the longitudinal fractional 
momenta x = (xi,x 2 , ....x n ) and virtualness up to Q 2 within the hadron and, in general, 
includes all Fock states with quantum numbers of the hadron. But only the lowest Fock 
state (gig2-for mesons, uud-for proton, etc.) contributes to the leading scaling behavior, 
other Fock state contributions are suppressed by powers of 1/Q 2 - In our work, we have 
restricted ourselves to considering the lowest Fock state for a meson. Then x = X\,x 2 
and Xi + x 2 = 1. This approach can be applied not only to the investigation of exclusive 
processes, but also to the calculation of higher twist corrections to some inclusive processes. 
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The qiq 2 spin state used in computing Tjj may be written in the form 



E 



u s1 {xiPm)v S2 (x2Pm) 




^f, p L helicityO, 
^T £ -^f, Pt helicity±l, 



V2 ' 



(2.2) 



where e± = t(1/v^)(0, 1, ±i, 0) in a frame with (pm)i,2 — and the N* project out 
a state of spins s, and pu is the four-momentum of the final meson. In our calculation, 
we have neglected the pion and the proton masses. Turning to extracting the contributions 
of the high twist subprocesses, there are many kinds of leading twist subprocesses in pp 
collisions as the background of the high twist subprocess qi + q 2 — > n + (or n~) + 7, such as 

q + q -> 7 + 9(9 -»■ tt + (*■"))» 9 + 5 -> 7 + 5(5 -> 7T + (7r~))» 9 + 5 -> 7 + 50(5 -> vr + (7r-)) and 
so on. The contributions from these leading twist subprocesses strongly depend on some 
phenomenological factors, for example, quark and gluon distribution functions in proton 
and fragmentation functions of various constituents etc. Most of these factors have not 
been well determined, neither theoretically nor experimentally. Thus they cause very large 
uncertainty in the computation of the cross section of process pp — > n + (or n~) + 7 + X. In 
general, the magnitude of this uncertainty is much larger than the sum of all the high twist 
contributions, so it is very difficult to extract the high twist contributions. 

The Mandelstam invariant variables for subprocesses q\ + q\ — > 7r + (7r~) + 7 are defined 

as 



In our calculation, we have also neglected the quark masses. We have aimed to calculate 
the pion production cross section and to fix the differences due to the use of various pion 
model functions. We have used five different functions: the asymptotic wave function ASY, 
the Chernyak-Zhitnitsky [2,5], ^2,-^3 model functions [7,8] and the wave function in which 
two non-trivial Gegenbauer coefficients a 2 and a 4 have been extracted from the CLEO data 
on the 77* — > 7T° transition form factor [20]. In ref. [20], the authors have used the QCD 
light-cone sum rules approach and have included into their analysis the NLO perturbative 
and twist-four corrections. They found that in the model with two nonasymptotic terms, at 
the scale pL = 2AGeV. 



s = 



(P1+P2) 2 , i=(p 1 -p 7T ) 2 , u = ( Pl - Py ) 2 . 



(2.3) 



$asy(x) = VSf n x(l - X), <S> C z(x, fl%) = 5^ asy (2x - if, 
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<5>p 2 {x,hI) = ®asy{x) [-0.1821 + 5.91(2x - l) 2 ], 
$P 3 (a;, / u 2 ) = $ asj/ (a;)[0.6016 -4.659(2a;- l) 2 + 15.52(2x - l) 4 ], 

$cW^o) = W^t 1 + 0.19C 2 3/2 (2x - 1) - 0.14C 4 3/2 (2:r - 1)], (2.4) 

where /^=0.0923 GeV is the pion decay constant. Here, we have denoted by x = x±, the 
longitudinal fractional momentum carried by the quark within the meson. Then, X2 = 1 — x 
and x\ — x 2 — 2x — 1. The pion wave function is symmetric under replacement x\ — x 2 <-> 
x 2 — x\. The values of the pion wave function moments < £ n > are defined as 

< C>= f (2.5) 

Here, is the model function without f n and £ = Xi — x 2 . The pion wave function 

moments have been calculated by means of the QCD sum rules method by Chernyak and 
Zhitnitsky at the normalization point /io = 0.5GeV. They are equal to 

< e >, = 1, < e >, = 0.44, < £ 4 >,= 0.27 (2.6) 
The Chernyak-Zhitnitsky pion model wave function has the following moments 

< £° > w = 1, < e 2 >, = 0.43, < > w = 0.24 (2.7) 

It is interesting to note that the corresponding moments of the asymptotic wave function 
differ considerably from those in Eqs.(2.6), (2.7) 

< £° > w = 1, < e >,= 0.20, < e >, = 0.086 (2.8) 

This means that the realistic pion wave function is much wider than the asymptotic one 
[5,21]. We have also used two other functions P 2 and P 3 . The wave function P 2 is a 
quadratic polynomial. Its free parameter, which is defined as the ratio of the coefficients of 
the constant and quadratic terms, is fixed to give the best fit to the values of the 2nd and 
4th moments of the wave function. The model function P 2 has the same shape as the CZ 
one, the difference being in the constant term. The function P 3 is a polynomial in which 
the number of free parameters is equal to the number of independent moments. These 
parameters are completely fixed by the same moment of the pion wave function. P 3 has a 
form different from that of the CZ and P 2 ones. All these wave functions, of course, have 
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the same zeroth moment by construction. Thus, the comparison of the predictions obtained 
by using the CZ and P2 ones with the ones obtained by means of P3 enables us to determine 
the sensitivity of the predictions to the form of the wave function. 
The model functions can be written as 

®asy(x) = V3f w x(l ~X), 



C 3 /2 (2x-l) + ^C? z (2x-l) 



,3/2, 



$p 2 (x,^) = $ asy (x) [c 3/2 (2x - 1) + 0.788C 2 3/2 (2x - 1) 
®p 3 (x,i4) = $ asy (x) [c 3/2 (2x - 1) + 0.7584C 3/2 (2a; - 1) + 0.3942C 3/2 (2a; - 1)] , 
®CLEo(x,fil) = ^asy(x) [c 3/2 (2x - 1) + 0.19C 3/2 (2:r - 1) - O.UC 3 4 /2 (2x - 1)] , 



(2.9) 



-y3/2 
) 

,3/2 



( 'o /J (2.r - 1 ) = 1, C 3/2 (2x - 1) = -(5(2* - l) 2 - 1), 
15 



Ct'\2x - 1) = ^-(2l(2x - l) 4 - 14(2x - l) 2 + 1). 
8 

It may be seen that the pion wave function extracted from the experimental data depends on 
the methods used and their accuracy. Although one may claim that the meson wave function 
is a process-independent quantity, describing the internal structure of the meson itself, the 
exploration of different exclusive processes with the same meson leads to a variety of wave 
functions. This means that the methods employed have shortcomings or do not encompass 
all the mechanisms important for a given process. Such a situation is pronounced in the 
case of the pion. It is known that the pion wave function (distribution amplitude- DA) can 
be expanded over the eigenfunctions of the one-loop Brodsky-Lepage equation, i.e., in terms 
of the Gegenbauer polynomials {C 3 J 2 (2x — 1)}, 



$ir(x,Q 2 ) = $ asy (x) 



1+ £ a n (Q 2 )C 3 J 2 (2x-l) 



n=2A. 



(2.10) 



The evolution of the wave function (DA) on the factorization scale Q 2 is governed by the 
functions a n (Q 2 ), 



/,o2\ / 2\ a s(Q 2 ) 

a n {Q ) = a„(/i ) — j—^r 
_<Wo) 

72 50 74 364 



(2.11) 



In Eq.(2.11), {7„} are anomalous dimensions defined by the expression, 



7n = Cf 



1 



n+1 -, 



(n + l)(n + 2) ~j 



(2.12) 



The constants a n (/ig) = a° are input parameters that form the shape of the wave functions 
and which can be extracted from experimental data or obtained from the nonperturbative 
QCD computations at the normalization point The QCD coupling constant a s (Q 2 ) at 
the two-loop approximation is given by the expression 



a s (Q 2 ) = 



An 



(3 ln{Q 2 /K 2 



2p 1 lnln(Q 2 /A 2 ^ 



(2.13) 



Here, A is the QCD scale parameter, /3q and f3\ are the QCD beta function one- and two-loop 
coefficients, respectively, 

2 19 
Po = 11 - -n/, A = 51 - —n f . 

The cross section for the high twist subprocess is given by the expression 



da,, . ^ 8ir 2 a E C F [D(t,u)\ 

— Is.t.u) = — — - 

dV ! 27 s 3 



1 1 

(j 2 + ¥ 



(2-14) 



where 



D(i, it) = e\t l dx\ 



a s (Q 2 )^uQj) 
1 — X\ 



+ 6 2 U / dX\ 
'0 



1 — X\ 



(2.15) 



where Q\ = {x\ — l)u, Q\ = —xii, represents the momentum squared carried by the hard 
gluon in Fig.l, ei(e 2 ) is the charge of qi(q 2 ) and Cf = §■ The high twist contribution to the 
large-^T pion production cross section in the process pp — > 7r + (7r~) + 7 is [22-24] 

j:Z T = E^-= [ [ dx 1 dx 2 G qi/hl (x 1 )G q2/h2 (x 2 )-^(qq^ 7:^)5(8 + t + u). (2.16) 
" V Jo Jo ft dt 



irE 



da 



da 

d 3 p dydp\ ' 

S = X±X 2 S, 



t = x±t, 



u = x 2 u, 
t = —m T y/se~ y = — pT\fse~ v , 



(2-17) 
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Xi 



x 2 



u = —rriT^/se y = —pT\/~se y ) 
x 2 u x 2 p T y/se y 



x 2 s + t x 2 s - prVse-y' 

X\t X\PT\fs e ^ V 



XiS + U X\S — PT^fsG y ' ! 

where mr - is the transverse mass of pion, which is given by 

2 2,2 

m T = m + p T . 

For a full discussion, we consider a difference A HT between the high twist cross section 
combinations H H T and Y, HT 

A HT = Sf+ T _ sf _T = E^{pp - 7T+7) - ^-^(PP - 7T-7). (2.18) 

We have extracted the following high twist subprocesses contributing to the two covariant 
cross sections in Eq.(2.16) 

da 1 - , da 2 - , da 3 _ da 4 _ 

—£r{ud -> vr + 7), —(d u ->7r"7), -^-(ud ^ ir 7), -^(du ^ n 7), (2.19) 

By charge conjugation invariance, we have 

da 1 . — , . da 3 . _ . da 2 . — , dcr 4 , _ , . 

— — [ud — > 7r 7) = — — (im — > Ti 7), ana — — (au — > 7r 7) = — — (au — > 7r 7). (2.20) 
dt dt dt dt 

III. CONTRIBUTION OF THE LEADING TWIST DIAGRAMS 

Regarding the high twist corrections to the pion production cross section, a comparison 
of our results with leading twist contributions is crucial. The leading twist subprocesses 
for the pion production are quark-antiquark annihilation qq — > (77, g — > 7i + (n~), shown in 
Fig. 2. The corresponding cross section is easily verified as[22] 

- i {qq^gq) = -na E a s {Q)-[ J + i y (3.1) 

For the leading-twist contribution, we find 

d(j s dcr 

Sm = E ^p~ = Yl / dx i dx 2dzG qi/hl (x 1 )G q2/h2 (x 2 )D n g (z) — ^(qq^g 1 )5(s + i + u), 

q 

(3.2) 

where 

s = Xix 2 s, t = — , u = , z = . (3.3) 

Z Z X\X 2 S 



Dg(z) = Dg (z) = Dg (z) represents the gluon fragmentation function into a meson con- 
taining a gluon of the same flavor. In the leading twist subprocess, n meson is indirectly 
emitted from the gluon with fractional momentum z. The S function may be expressed in 
terms of the parton kinematic variables, and the z integration may then be done. The final 
form for the cross section is 

^ L m = E ~p- = / dxi / dx 2 G qi/hl (x 1 )G q2/h2 (x 2 )Dg(z)x 



KZdV ' ^J XJ . J X2 . -(X 1 t + X 2 u) 7T dt K 



^lmin x 2n 

(3.4) 



IV. NUMERICAL RESULTS AND DISCUSSION 



In this section, the numerical results for higher twist effects on the dependence of the 
chosen meson wave functions in the process pp — > n + (or are discussed. We have 

calculated the dependence on the pion wave functions for the high twist contribution to the 
large-pr pion production cross section in the proton-proton collision . In the calculations, 
the asymptotic & asy , Chernyak-Zhitnitsky §cz, two other pion model functions, <3>p 2 , $p 3 
and also, the pion wave function, from which two non-trivial Gegenbauer coefficients a 2 
and a 4 have been extracted from the CLEO data on the 7r°7 transition form factor have 
been used[20]. In the ref.[20], authors have used the QCD light-cone sum rules approach 
and included into their analysis the NLO perturbative and twist-four corrections. For the 
high twist subprocess, we take qi + q 2 (Q1Q2) + 7 and we have extracted the following 
four high twist subprocess ud — > du — > 7T + 7, ud — > 7r~7, du — > 7r~7 contributing to 

pp — > 7r + (or 7r~)7 cross sections. For the dominant leading twist subprocess for the pion 
production, we take the quark- ant iquark annihilation qq — > (77, in which the n meson is 
indirectly emitted from the gluon. As an example for the quark distribution function inside 
the proton, the MRST2003c package [25] has been used. The gluon fragmentation function 
has been taken from [26]. The other problems dealth with are the choice of the QCD scale 
parameter A and the number of the active quark flavors rif. The high twist subprocesses 
probe the meson wave functions over a large range of Q 2 squared momentum transfer, carried 
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by the gluon. Therefore, phenomenologically, in the given diagram in Fig 1, for the center- 
of-mass energy y/s = 63GeV, we take Q 2 = p\. However, for the center-of-mass energy 
yfs = 630Ge^, we take Q 2 = {x\ — l)u, Q\ = —X\t , which we have obtained directly 
from the high twist subprocesses diagrams. The same Q 2 has been used as an argument 
of a s (Q 2 ) in the calculation of each diagram. The results of our numerical calculations are 
plotted in Figs. 3-18. Figs. 3-5 show the dependence of the differential cross sections of the 
high twist Sf7, leading twist plus high twist S^=S^+SfJ and ratio R = EfJ/S^J as 
a function of the pion transverse momentum px for five different meson wave functions. 
As shown in Figs. 3-4, the high twist differential cross section is monotonically decreasing 
with an increase in the transverse momentum of the pion. As seen from Figs. 3-4, in all 
wave functions of the mesons, the dependencies of the high twist cross sections on the pt 
transverse momentum of the pion demonstrate the same behavior. On the other hand, 
the higher twist corrections are very sensitive to the choice of the pion wave function. We 
should note that the magnitude of the high twist cross section in the pion wave function 
&cleo( x ,Q 2 ) case is very close to the asymptotic wave function <& aS y(x) case. In Fig. 5, the 
ratio R = S^j r /E^+ is plotted at y = as a function of the pion transverse momentum p? 
for the different pion wave functions. First of all, it is seen that the values of R for fixed 
y and depend on the choice of the pion wave function. Also, the distinction between 
R($ asy (x)) with R($cleo(x,Q 2 )), R($cz(x,Q 2 )), R($ P2 (x,Q 2 )) and R($p 3 (x,Q 2 )) have 
been calculated. We have found that the distinction R($ asy (x)) and R(&cleo( x ,Q 2 )) is 
small, whereas a distinction between R($ asy (x)) with R(§cz( x ,Q 2 )), R{^p 2 (. x iQ 2 )) an d 
R(§p 3 (x, Q 2 )) is significant. For example, in the case of \fs = 63GeV, y = 0, the distinction 
between R(<$> asy (x)) with R{<S>i{x,Q 2 )) (i = CLEO,CZ, P 2 , P 3 ) is shown in Tabled 

Thus, the distinction between R($ asy (x)) and R(§cleo(x, Q 2 )) is maximum at pr = 
6GeV/c, but the distinction between R($ asy (x)) with R(§ C z(x, Q 2 )), R(®p 2 ( x ,Q 2 ))i 
R(&p 3 (x,Q 2 )) is maximum at px = 2GeV/c and decreases with an increase in p T . Such 
a behavior of R may be explained by reducing all moments of the pion model wave func- 
tions to those of $ asj/ (x) for high Q 2 . In Fig. 6, we show the difference of the A^ T =S^ r -S^ T 
high twist cross section as a function of the pion transverse momentum p? for five different 
pion wave functions. As seen from Fig.6, when the transverse momentum of the pion is in- 
creasing, the difference A^ T high twist cross section is monotonically decreasing. As shown 
in Fig.6, the difference A^ T high twist cross section for pion wave function &cleo( x , Q 2 ) 
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is very close to the asymptotic wave function & asy (x). In Fig. 7, we show the difference 
A£ T =£^+-£^T leading twist and sum of differences leading and high twist A t ° t =A% r +A% T 
cross section for five pion wave functions as a function of the pion transverse momentum p?- 
As in Fig. 6, the A^ T and A 1 " 1 are monotonically decreasing when the transverse momentum 
Pt of the pion is increasing. But, the distinction between the difference A^ T leading and 
A t ° t =A^ T +A^ T sum of the difference leading and the high twist cross section is not evident. 
In Fig. 8, the ratio r = A^ T /A^ T is plotted at y = as a function of the pion transverse 
momentum p? for five pion wave functions. As shown in Fig.8, the values of r for fixed y 
and y/s depend on the choice of pion wave function as in Fig.5. Also, we have calculated 
the distinction between r($ asy (x)) with v(^cleo(x, Q 2 )), r($cz(x,Q 2 )), r($p 2 (2;,Q 2 )) and 
r($p 3 (x, Q 2 ))- For example, in the case of y/s = 63GeV, y = the distinction between 
r{<$> asy {x)) with r(^i(x,Q 2 )) (i=CLEO, CZ, P 2 , P 3 ) is presented in TableHH 

We have obtained very interesting results. The calculations show that the ratio 
Q 2 ))/ R(§ asy (x)), (i=CLEO, CZ, P 2 , P3) of all the transverse momentum pt of 
the pion is equal to the ratio r($i(x, Q 2 ) /r($ asy (x)), (i=CLEO, CZ, P 2 , P?)- In Figs.9-10, 
we have depicted £^+" and A^ T =S^ r -E^ r as a function of the rapidity y of the pion at 
y/s = 63GeV and pt = 5GeV/c. As we are now in the high energy region, the change of 
the rapidity to determine these relations is given by —ln(y/s/p T ) < y < /n(-y/i/pr)- At 
y/s = 63GeV and px = 5GeV/c, the pion rapidity lies in the region —2.52 < y < 2.52. 
As seen from Figs.9-10, in the region (—2.52 < y < 0.85), the high twist cross section for 
asymptotic $> asy (x) and §cleo{%,Q 2 ) wave functions increase and it has a maximum ap- 
proximately at one point y = 0.85. After that, the cross sections monotonically decrease 
with an increase in the y rapidity of the pion. As seen from Figs.9-10, the high twist ^+ 
and difference high twist A^ T cross sections are very sensitive to the choice of meson wave 
functions. Also, as shown in Figs.9-10, in the regions —2.52 < y < — 1 and 1 < y < 2.52, 
the high twist cross section for the $cleo(x, Q 2 ) is very close to the $ asy (x) wave func- 
tion case. We have also carried out comparative calculations in the center-of-mass energy 
y/s = 630GeV. Figs. 11-13 show the dependence of the differential cross sections of the 
high twist SfJ, leading twist plus high twist £^=£^+£^7 and ratio R = EfJ/E^ as 
a function of the pion transverse momentum p? for five different meson wave functions. 
As shown in Figs. 11-12, the high twist differential cross section is monotonically decreas- 
ing when the transverse momentum of the pion is increasing. As seen from Figs. 11- 12, 
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for all wave functions of mesons, the dependencies of the high twist cross sections on the 
transverse momentum of the pion demonstrate the same behavior. But, the higher twist 
corrections are very sensitive to the choice of the pion wave function. We should note that 
the magnitude of the high twist cross section for the pion wave function §cleo{ x , Q 2 ) is 
very close to the asymptotic wave function $ asy (:r) case. In Fig. 13, the ratio R = E^/E^ 
is plotted at y = as a function of the pion transverse momentum px for the different 
pion wave functions. First of all, it may be seen that the values of R for fixed y and ^/s 
depend on the choice of the pion wave function. As in center-of-mass energy y/s = 63GeV, 
we have also calculated the difference between R(§ asy (x)) and four other wave functions 
i.e. R($cleo(x,Q 2 )), R(®cz(x,Q 2 )), R($p 2 (x,Q 2 )), and R($ P3 (x,Q 2 )). We have found 
that the distinction between R(& asy (x)) and R{^cleo{ x ^Q 2 )) is ver Y small, whereas the 
distinction between R($ asy (x)) and R($ C z(x, Q 2 )), R(®p 2 (x, Q 2 )) and R($p 3 (x, Q 2 )) is sig- 
nificant. In order to demonstrate this, in the case of \fs = 630GeV, y = the distinction 
between R(<$> asy (x)) with R($i(x,Q 2 )) (i=CLEO, CZ, P 2 , P 3 ) is shown in Table ITTT1 

Thus, the distinction between R(Q asy (x)) and R(§cleo(%, Q 2 )) is maximum at pp = 
QOGeV/c, but the distinction between R($ asy (x)) with R(§ C z(x, Q 2 )), R(®p 2 ( x ,Q 2 ))i 
R(&p 3 (x,Q 2 )) is maximum at pt = 20GeV/c and decreases with an increase in pp. In 
Fig. 14, we have shown the difference A^ T =Y^+ -Y^I high twist cross section as a function 
of the pion transverse momentum pp, for five different pion wave functions. As seen from 
Fig. 14, as the transverse momentum of pion, pt, increases, the difference of the high twist 
cross section, A^ T , monotonically decreases. As shown in Fig. 14, A^ T for the $>cleo(%, Q 2 ) 
pion wave function is very close to the asymptotic wave function, § asy (x). In Fig. 15, we have 
shown A^ T =£^+-£^T i.e. the leading twist cross section, and A t ° t =A% r +A% T i.e. the sum 
of the leading and the high twist cross sections, for five pion wave functions, as a function of 
the pion transverse momentum, pp. As in Fig. 7, the A 7T and A^°* cross sections monotoni- 
cally decrease when the transverse momentum pt of the pion increases. But the distinction 
between A 7T and A^°*=A^ T +A^ T is not evident. In Fig.16, the ratio r = A^ T / A 7T is 
plotted at y = as a function of the pion transverse momentum, pt, for five pion wave 
functions. As shown in Fig.16, the values of r for fixed y and \fs depend on the choice of 
pion wave function similar to Fig.8. We have also calculated the distinction between r(Q asy ) 
and r($cLEo( x i Q 2 )> r (^cz(x,Q 2 )), r($p 2 (x,Q 2 ) and r(§ P . 3 (x, Q 2 ). As an example, in the 
case of y/s = 630GeV, y = the distinction between r(§ asy (x)) and r($ i (x, Q 2 )) (i=CLEO, 
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CZ, P 2 , P 3 ) is shown in Table HVl 

As seen from calculations with increasing center-of-mass energy from a/s = 63GeV to 
y/s = 630Ge^, the distinction between R and r decreases for all pion wave functions. In 
Figs. 17-18, we have depicted ^+ and A^ T =H^+-H^T as a function of the rapidity y of 
pion at the \/s = Q30GeV and pr = 50GeV/ c. Since we are now in the high energy region, 
the change of the rapidity to determine these relations is —ln{\fs/pT) < y < ln(\fs / 'pr) ■ At 
y/s = 630GeV and pt = 50GeV/c the pion rapidity also lies in the region —2.52 < y < 2.52. 
As seen from Figs. 17-18, in the region (—2.52 < y < 0.85), the high twist cross section for 
the asymptotic <3> asy (x) and $cleo(%, Q 2 ) wave functions increases and it has a maximum 
approximately at one point y = 0.85. After that, the cross sections monotonically decrease 
with an increase in the y rapidity of the pion. As in Figs. 9-10, in the regions —2.52 < y < — 1 
and 1 < y < 2.52, the high twist cross section for &cleo( x i Q 2 ) is ver Y close to the Q a sy(%) 
wave function. 



V. CONCLUDING REMARKS 



In this work, we have calculated the higher twist contribution to the large-p^ pion pro- 
duction cross section to show the dependence on the chosen meson wave functions in the 
process pp — > n + (or 7r~)7. In our calculations, we have used the asymptotic § a sy, Chernyak- 
Zhitnitsky $cz, two other pion model functions, <3>p 2 , $p 3 and also, the pion wave function, 
in which the coefficients a 2 and have been extracted from the CLEO data on the 7r°7 tran- 
sition form factor used. For the high twist subprocess, we have taken qi+q 2 (<?i<?2)+7- We 
have extracted the following four high twist subprocesses ud 7r + 7, du — > vr + 7, ud — > 7r~7, 
du — > 7r~7, contributing to pp — > n + (orn^)^ cross sections. As the dominant leading twist 
subprocess for the pion production, we have taken the quark-antiquark annihilation qq — > y-f, 
where the 7r meson is indirectly emitted from the gluon. The results of our numerical cal- 
culations have been plotted in Figs. 3-18. As shown in Figs. 3-4 and Figs. 11-12, the high 
twist differential cross section monotonically decrease when the transverse momentum of 
the pion increases. As seen from Figs. 3-4 and Figs. 11- 12 in all wave functions of mesons, the 
dependencies of the high twist cross sections on the px transverse momentum of the pion 
demonstrate the same behavior. But, the higher twist corrections are very sensitive to the 
choice of the pion wave function. It should be noted that the magnitude of the high twist 
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cross section for the pion wave function &cleo(%, Q 2 ) is very close to the asymptotic wave 
function $ asy (x). 

In Figs. 5 and 13, the ratio R = S^ r /E^+ has been plotted at y = as a function 
of the pion transverse momentum, p^, for the different pion wave functions. It may be 
observed that the values of R for fixed y and y/s depend on the choice of pion wave func- 
tion. Within this context, we have also calculated the distinction between R($ asy (x)) and 
R($cleo(x,Q 2 )), R($cz(x,Q 2 )), R($p 2 (x,Q 2 )) and R(<S> Pi (x,Q 2 )). We have ultimately 
found that the difference between R(& asy (x)) and R(&cleo( x ,Q 2 )) is small, whereas a dis- 
tinction between R(§ asy (x)) with R(& C z(x, Q 2 )), R(®p 2 ( x ,Q 2 )) an d -R($p 3 (x, Q 2 )) is sig- 
nificant. In Figs. 6 and 14, we have shown the difference high twist cross section, A^ T =S^ r_ 
S^ T , as a function of the pion transverse momentum, pt, for five different pion wave func- 
tions. As seen from Figs. 6 and 14 when the transverse momentum of the pion increases, the 
difference of the high twist cross section, A^ T , monotonically decreases. As shown in Figs. 6 
and 14, A^ T in the $cleo(x, Q 2 ) pion wave function is very close to that of the asymptotic 
wave function $ asy (x). 

In Figs. 8 and 16, the ratio r = A^ T /A^ T is plotted at y = as a function of the 
pion transverse momentum, pt, for five pion wave functions. As shown in Figs. 8 and 16, the 
values of r for fixed y and \/s depend on the choice of the pion wave function similar to Fig. 5 
and 13. We have also calculated the distinction between r(Q asy (x)) and r(Q C LEo(%, Q 2 )), 
r($cz(x,Q 2 )), r(<frp 2 (x,Q 2 )) and r($p 3 (:r, Q 2 )). For all transverse momentum p T of the 
pion in the center-of-mass energy \/s = 63GeV and also in 630GeV, we have obtained the 
following interesting relation: 

R { H^Q 2 ))^r(H X ,Q 2 )) {i = CLE o,CZ,P 2 .P 3 ) 

R($ asy (x)) r($ asy (x)) K ,,23) 

In Figs. 9- 10 and Figs. 17- 18, we have depicted high twist and difference high twist cross 
sections A^ T =E^ r -E^ T as a function of the rapidity y of the pion at y/s = 63GeV, Q30GeV 
and pt = 5GeV/ c, 50GeV/ c, respectively. As we are now in the high energy region, the 
change of the rapidity of this relation may be expressed as folows: —ln(y/s/pT) < y < 
ln(y/s/p T ). At v/i = 63GeV, p T = hGeV/c and yfs = 630GeV, p T = 50GeV/c the pion 
rapidity lies in the region —2.52 < y < 2.52. As seen from Figs. 9-10 and Figs. 17-18 in 
the region (—2.52 < y < 0.85), the high twist cross section for the asymptotic & aS y(x) and 
&cleo(x,Q 2 ) wave functions increases and it has a maximum approximately at one point 
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y = 0.85. After that, the cross sections monotonically decrease with an increase in the y 
rapidity of the pion. As seen from Figs. 9-10 and Figs. 17-18, the and the A^ T cross 
sections are very sensitive to the choice of the meson wave functions. Also, as shown in 
Figs. 9, 10, 17 and 18 in the regions —2.52 < y < — 1 and 1 < y < 2.52 the high twist cross 
section for the §cleo{ x , Q 2 ) is ver y close to the asymptotic wave function $> asy (x) case. Our 
investigation enables us to conclude that the high twist pion production cross section in the 
proton-proton collisions depends on the form of the pion model wave functions and may be 
used for their study. Further investigations are needed in order to clarify the role of high 
twist effects in QCD. 

Acknowledgments 

Two of authors, A. I. Ahmadov and I. Boztosun are grateful to TUBITAK Grant-2221 
(BAYG) as well as TUBITAK Grant: TBAG-2398. One of us, A. I. Ahmadov is also grateful 
to NATO Reintegration Grant-980779 



VI. REFERENCES 



A.V. Radyushkin, Dubna preprint P2-10717, 1977, hep-ph/0410276. 



V. L. Chernyak and A. R. Zhitnitsky, Nucl. Phys. B201, 492 (1982); V. L. Chernyak, A. 
R. Zhitnitsky and I. R. Zhitnitsky, Nucl. Phys. B204, 477 (1982). 
V. L. Chernyak and A. R. Zhitnitsky, Nucl. Phys. B246, 52 (1984). 
I. D. King, C. T. Sachrajda, Nucl. Phys. B279, 785 (1987). 
V. L. Chernyak, A. R. Zhitnitsky, Phys. Rept. 112, 173 (1984). 

S. V. Mikhailov and A. V. Radyushkin, JETP Lett. 43, 712 (1986); Sov. J. Nucl. Phys.49, 

494 (1989); Phys. Rev. D45, 1754 (1992). 

V. M. Braun and I. E. Filyanov, Z. Phys. C44, 157 (1989). 

G. R. Farrar, K. Huleihel and H. Zhang, Nucl. Phys. B349, 655 (1991). 

V. L. Chernyak, A. A. Ogloblin and I. R. Zhitnitsky, Z. Phys. C42, 583 (1989). 

G. R. Farrar, H. Zhang, A. A. Ogloblin and I. R. Zhitnitsky, Nucl. Phys. B311, 585 (1988). 

M. Anselmino, P. Kroll and B. Pire, Z. Phys. C36, 89 (1987). 



16 



[12] A. V. Radyushkin and R. Ruskov, Phys. Lett. B374, 173 (1996); Nucl. Phys. B481, 625 
(1996). 

The CLEO Collaboration (J.Gronberg et.al), Phys. Rev. D57, 33 (1998). 
V. Yu. Petrov, M. V. Polyakov, R. Ruskov, C. Weiss and K. Goeke, Phys. Rev. D59, 114018 



(1999); hep-ph/9807229. 

J. Qiu and G. Sterman, Nucl. Phys. B353, 105, 137 (1991). 

E. L. Berger and S. J. Brodsky, Phys. Rev. Lett. 42, 940 (1979); E. L. Berger, Z. Phys. C4, 
289 (1980). 

R. Doria, J. Frenkel and J. C. Taylor, Nucl. Phys. B168, 93 (1980); F. T. Brandt, J. Frenkel 
and J. C. Taylor, Nucl. Phys. B312, 589 (1989). 

K. Kastella, G. Sterman and J. Milana, Phys. Rev. D39, 2586 (1989). 
G. P. Lepage and S. J. Brodsky, Phys. Rev. D22, 2157 (1980). 
A. Schmedding and O. Yakovlev, Phys. Rev. D62, 116002 (2000), |hep^ph/9905392 
V. L. Chernyak, Preprint TPI-MINN-91-47-T, December 1991. 
J. F. Owens, Rev. Mod. Phys. 59, 465 (1987). 

R. Kh. Muradov and A. I. Ahmadov, Central. Eur. J. Phys. 3(3) (2005) 433, |hep-ph/0412227[ 
A. I. Ahmadov, I. Boztosun, R. Kh. Muradov, A. Soylu and E. A. Dadashov, hep-ph/0604172, 
MRST2003c.f can be obtained from http://durpdg.dur.ac.uk/hepdata/pdf.html, See also, 
A. D. Martin, R. G. Roberts, W. J. Stirling and R. S.Thorne, |hej>ph /0307262 ; R. S.Thorne, 
|hep-ph/0309343[ 

[26] B. A. Kniehl, G. Kramer, B. Potter, Nucl.Phys. B582, 514 (2000), |hep^ph/0010 289 



17 



PT ,GeV/c 


R(®cleo(x,Q 2 )) 


R(®cz(x,Q 2 )) 


R(4>P 2 (x,Q 2 )) 


R(4>P 3 (x,Q 2 )) 


nyVasyiyX)) 


n^asy(X)) 


H(^a3y(X)) 


ti[^asy\X)) 


2 


1.496 


8.435 


10.569 


16.525 


6 


3.358 


3.888 


4.617 


3.001 


20 


0.9476 


0.382987 


0.302073 


0.133824 



TABLE I: The distinction between R(<f> asy (x)) with R(<S>i(x,Q 2 )) (i=CLEO, CZ, P 2 , P 3 ) at cm. 
energy \fs = QSGeV . 



p T , GeV/c 


r($CLEo(x,Q 2 )) 


r\3>Cz(x,Q 2 )) 


r(<3>p 2 (x,Q 2 )) 


r(4>p 3 (x,Q 2 )) 


r(® asy (x)) 


r(® asy (x)) 


r($ asy (x)) 


r(*««v(*)) 


2 


1.496 


8.435 


10.568 


16.524 


6 


3.358 


3.887 


4.616 


3.001 


20 


0.95 


0.38315 


0.302109 


0.13384 



TABLE II: The distinction between r(<S> asy (x)) with r(<S>i(x,Q 2 )) (i=CLEO, CZ, P 2 , P 3 ) at cm. 
energy \fs = 63Gey. 



Pt, GeV/c 


R($cleo(x,Q 2 )) 

R($a S y(x)) 


R($cz(x,Q 2 )) 

R($a S y(x)) 


R($p 2 (x,Q 2 )) 

R(®asy(x)) 


R(4>P 3 (x,Q 2 )) 

R{$asy(x)) 


20 


1.531 


5.582 


6.815 


9.63 


60 


2.202 


2.591 


2.961 


2.262 


200 


0.926 


0.605 


0.544 


0.43 



TABLE III: The distinction between R(<$> asy (x)) with R(<S>i(x,Q 2 )) (i=CLEO, CZ, P 2 , P 3 ) at cm. 
energy y/s = 630GeV. 
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PT ,GeV/c 


r(<S>CLEo(x,Q 2 )) 


r(^cz(x,Q 2 )) 


r($P 2 (x,Q 2 )) 


r(4>p 3 (x,Q 2 )) 






r(1?asy(X)) 




20 


1.531 


5.583 


6.817 


9.632 


60 


2.202 


2.591 


2.961 


2.262 


200 


0.926 


0.605 


0.544 


0.43 



TABLE IV: The distinction between r(<$> asy (x)) with r(<S>i(x,Q 2 )) (i=CLEO, CZ, P 2 , P 3 )) at cm. 
energy y/s = 630GeV. 
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FIG. 4: The sum of the leading and the high twist ir + production cross sections S*°* = S L J + S 
as a function of the pt transverse momentum of the pion , at the cm. energy yfs = 63GeV . 
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FIG. 5: Ratio i? = S^ r /S^+ , where the leading and the high twist contributions are calculated 
for the pion rapidity y = at the cm. energy ^fs = 63GeV ', as a function of the pion transverse 
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FIG. 6: The difference of the high twist cross section, A^ T =S^ r -S^ T , as a function of the pion 
transverse momentum, pj>, at the cm. energy yfs = &?>GeV . 
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FIG. 7: The sum of the difference leading and high twist cross sections A t ° t =A^ T +A^ T as a 
function of the pion transverse momentum px at the cm. energy yfs = 63GeV. 
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FIG. 8: The ratio r = A^ T / A^ T , where the leading and the high twist contributions are calculated 
for the pion rapidity y = 0, at the cm. energy y/s = 63GeV, as a function of the pion transverse 
momentum pr- 
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FIG. 9: High twist ir + production cross sections as a function of the y rapidity of the pion at the 

transverse momentum of the pion pt = 5GeV, at the c.m.energy y/s = 63GeV. 
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FIG. 10: The difference of the high twist cross section, A^ T =S^ r -S^ T , as a function of the y 
rapidity of the pion at the transverse momentum of the pion pt = 5GeV, at the cm. energy 
Js = Q3GeV. 
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High twist 7r + production cross sections as a function of the pr transverse momentum of 
at the cm. energy y/s = 630GeV. 
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FIG. 12: The sum of the leading and the high twist 7r + production cross sections = +£^ r 
as a function of the px transverse momentum of the pion, at the cm. energy y/s = 630Gey. 
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FIG. 13: The ratio R = E^/E^, where the leading and the high twist contributions are cal- 
culated for the pion rapidity y = at the cm. energy *Js = 630Gey as a function of the pion 
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FIG. 14: The difference of the high twist cross section, A^ T =E^ r -E^ r , as a function of the pion 
transverse momentum, pr, at the cm. energy y/s = 630GeF. 
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FIG. 15: The sum of the difference of the leading and the high twist cross sections, 

A t ° t =A^ T +A^ T , as a function of the pion transverse momentum, pr, at the cm. energy 

y/s = 630GeV. 
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FIG. 16: The ratio r = A^ T / A^ T , where the difference of the leading and the high twist contri- 
butions are calculated for the pion rapidity, y = 0, at the cm. energy yfs = 630GeV, as a function 
of the pion transverse momentum p?. 
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FIG. 17: The high twist 7r + production cross sections as a function of the y rapidity of the pion at 
the transverse momentum of the pion px = 50GeF, at cm. energy y/s = 630GeV '. 
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FIG. 18: The difference of the high twist cross section, A^ T =S^ r -S^ r , as a function of the 
y rapidity of the pion at the transverse momentum of the pion px = 50GeV, at cm. energy 
Js = 630GeF. 
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